A Novel Route to Calculate the Length Scale for the Glass Transition in Polymers 
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The occurrence of glass transition is believed to be associated to cooperative motion with a 
growing length scale with decreasing temperature. We provide a novel route to calculate the size 
of cooperatively rearranging regions CRR of glass-forming polymers combining the Adam-Gibbs 
theory of the glass transition with the self-concentration concept. To do so we explore the dynamics 
of glass-forming polymers in different environments. The material specific parameter a connecting 
the size of the CRR to the configurational entropy is obtained in this way. Thereby, the size of CRR 
can be precisely quantified in absolute values. This size results to be in the range 1-^3 nm at the 
glass transition temperature depending on the glass-forming polymer. 

PACS numbers: 64.70.Pf, 83.80.Tc, 83.80.Rs, 77.22. Gm 



The nature of the glass transition is one of the 
most important unsolved problems in condensed mat- 
ter physics and research in this field has enormously 
intensified in the last decades due to its strong funda- 
mental as well as applicative implications. Among 
the peculiar phenomena displayed by glass-forming 
liquids, the super- Arrhenius temperature dependence 
of the viscosity and the structural correlation time is 
certainly one of the most intriguing. In this frame- 
work, more than forty years ago Adam and Gibbs 
[H theorized that such a pronounced temperature de- 
pendence of the structural correlation time is due to 
a cooperative process involving several basic struc- 
tural units forming cooperatively rearranging regions 
(CRR), which size increases with decreasing tempera- 
ture. Since then a great deal of theoretical approaches 
0, [3] as well as simulation studies [H and very re- 
cently experimental studies employing multipoint dy- 
namical susceptibilities @ have been devoted in the 
search of good candidates for CRR as well as its size 
and temperature dependence. All of these studies sug- 
gest that a growing correlation length with decreasing 
temperature of the order of several nanometers ex- 
ists. According to the Adam-Gibbs (AG) theory of the 
glass transition, the increase of the structural relax- 
ation time with decreasing temperature, accompanied 
by the growth of the cooperative length scale, is due to 
the decrease of the number of configurations the glass- 
former can access, namely the configurational entropy 
(S c ) of the system. The connection between the struc- 
tural relaxation time and S c is expressed by 



T exp[C/(TS c )], 



(1) 



where C is a glass-former specific temperature inde- 
pendent parameter and tq is the pre-exponential fac- 
tor. The dynamics of a large number of glass forming 
systems have been successfully described through the 
AG equation in both experiments for low molecular 
weight glass formers [f| and polymers Q, as well as 
simulations Apart from the relation between the 
relaxation time and the configurational entropy, the 
AG theory provides a connection between the number 
of basic structural units belonging to the CRR and 
the configurational entropy: N w S7 1 . Several re- 



cent simulations studies, where the cooperative length 
scale and the configurational entropy were obtained 
within the framework of the string-like motion and 
the potential energy landscape of the glass-former re- 
spectively, successfully tested the validity of this corre- 
lation Q . Being the number of particles proportional 
to the volume of the CRR, the characteristic length 
scale £ can be related to the configurational entropy 
by: 



C/2 



'1/3 



(2) 



where r c is the equivalent radius of the CRR and 
a is a material specific parameter. Here, the ex- 
ponent — 1/3 implies that the fractal dimension for 
the growth of CRR equals 3. This value means that 
CRR possess a compact shape growing in three direc- 
tions. This hypothesis has been rationalized within 
the framework of the random first order transition 
theory [Toj ] . Very recent molecular dynamics simu- 
lations in a Lennard-Jones liquid, performed taking 
into account the motion of all particles [ll[ , also sug- 
gest that the structural rearrangement involves the 
motion of compact regions. Despite the ability to de- 
scribe the temperature dependence of the structural 
relaxation and predict an increasing length scale with 
decreasing temperature, the glass-former specific pa- 
rameter a is not reliably obtained from the AG theory 
as highlighted by several authors 3, 12]. Hence the 
size of the CRR cannot be estimated a priori from the 
AG theory alone. In this work we apply the AG the- 
ory to polymer blends and polymer-mixture systems 
in order to provide a new route to precisely determine 
the parameter a and hence the size of CRR. To do 
this, we incorporate the concept of self-concentration 
in the AG theory, first proposed by Chung and Korn- 
field [l]| and later developed by Lodge and McLeish 
to describe the segmental dynamics of miscible 
polymer blends. This approach has been recently ex- 
ploited by us to describe the dynamics of miscible 
polymer blends and concentrated polymer solutions 
and provided an accurate description of the compo- 
nent segmental dynamics (the a relaxation) of these 
systems [lj| . The essential features of the concept of 
self-concentration are explained in Fig. 1 and can be 
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FIG. 1: (Color online) Schematic illustration of the self- 
concentration concept. 

summarized as follows: when a volume is centered on 
the basic structural unit of one of the polymers of a 
given blend, the effective concentration ((f> e g) will be 
in general different from the macroscopic one. Two 
extreme cases are possible: i) if this volume equals 
the volume of the basic structural unit, then c ff = 1; 
ii) if, on the other hand, this volume is large enough, 
e ff will be equal to the macroscopic concentration, 
4> e s = (/)■ The effective concentration is related to the 
self-concentration (<p s ) and to the macroscopic concen- 
tration through [l4|: </> c ff = <j) s + (1 — 4> s )4>- The self- 
concentration, namely the volume fraction occupied 
by the isolated chain containing the unit of reference, 
can be related to the radius of the volume through 
simple geometric considerations involving the Kuhn 
and the packing lengths (l6| . ft is worth remarking 
that, as the volume of CRR is expected to be neither 
large enough to make the concentration equal to the 
macroscopic one nor equal to that of the monomer 
{<f> < (peg < 1) [ii, 15, 1(1, the component dynamics 
will be intermediate between that of the pure poly- 
mer and the average structural dynamics of the mix- 
ture. This means that the self-concentration concept 
constitutes an extremely sensitive tool when exploring 
length scales of the order of those expected for CRR. 

In order to apply the AG theory to mixtures, S c , 
C, tq in Eqs. (1) and (2) have to be evaluated in the 
cooperative volume. This means that these quanti- 
ties are function of the effective concentration of the 
polymer under examination. In our previous works, 
we found that a linear combination of the parame- 
ters of the pure components of the mixture through 
the effective concentration allows a good description 
of dynamics data (IBj : 



X = <j) cS AXA + (1 - </>cffA)-X"E 



(3) 



Where X is either S c , C or logr . Eqs. (1), (3) 
and (4) can be fitted to experimental data having 
a as the only unknown parameter. As this is a 
glass-former specific parameter, the presence of the 
other component in CRR can in principle affect the 
value of a. This has been found to be the case es- 
pecially for polymer-solvent systems [Hj]. We have 



overcome this problem by studying the dynamics of 
polymer-polymer and polymer-solvent mixtures at el- 
evated concentration of the polymer for which a is 
being calculated. This allows a straightforward ex- 
trapolation of the a parameter to the pure polymer. 
Furthermore, it is noteworthy that for concentrated 
polymer mixtures both concentrations fluctuations, 
which might affect the dynamics, and deviations from 
ideality, which would induce the failure of Eq. (3), 
are minimized. We have thus obtained values of a 
for polystyrene (PS), poly(vinyl acetate) (PVAc) and 
used those already obtained for poly(o-chlorostyrene) 
(PoCIS) and polyvinyl methyl ether) (PVME) 
([22]). In such a way we have investigated four typi- 
cal glass-forming polymers with different rigidity. We 
have found that a single polymer specific value of a 
is obtained when extrapolating to 100% concentra- 
tion independently on the environment surrounding 
the polymer under consideration. The achievement 
of a temperature independent parameter results, via 
Eq. (3), in an increasing length scale for structural re- 
laxation with decreasing temperature of the order of 
nanometers. In particular, the resulting size of CRR 
was between 1 and 3 nm depending on the charac- 
teristic of the polymer. For a given polymer in the 
accessible temperature range the diameter of CRR 
only varies of about 20%. PVAc, poly(2,6-dimethyl- 
1,4-phenylene oxide) (PPhO), ortho-terphenyl (OTP) 
and toluene were purchased by Sigma- Aldrich. The 
molecular weight of the two polymers were respec- 
tively Af w = 83000g/mo; (M w /M m = 4.2) and M n = 
23000g/mol (M w /M n = 2.3). PS was purchased 
from Polymer Source Inc. The molecular weight was 
M n = 704005 jmol (M w /M n = 1.04). Bis-phenol- 
C-dimethylether (BCDE) was synthesized as describe 
elsewhere [17| . Highly concentrated PS and PVAc in 
different environments were investigated by means of 
broadband dielectric spectroscopy (BDS). PVAc and 
PS dynamics were investigated respectively in toluene 
and BCDE, and OTP and PPhO. Dielectric measure- 
ments were carried out using a high precision dielectric 
analyzer (ALPHA, Novocontrol GmbH) and a Novo- 
control Quatro cryosystem for temperature control 
with a precision of 0.1 K. Measurements were per- 
formed over a wide frequency (10~ 2 — 10 6 Hz) and 
temperature range in isothermal steps starting for the 
highest temperature. As a general rule the relevant re- 
laxation time at any temperature was obtained from 
dielectric relaxation spectra as the reciprocal of the 
angular frequency at the maximum of the permittiv- 
ity loss. Structural relaxation data for toluene [HI], 
BCDE 0, OTP and PPhO [H[ were taken from 
the literature. Calorimetric measurements were car- 
ried out on pure PVAc, PS and BCDE by means 
of the differential scanning calorimeter (DSC-Q1000) 
from TA-Instruments. Measurements were performed 
in temperature modulated mode with a average heat- 
ing rate of 0.1 K/min and amplitude of 0.3 K. Differ- 
ent oscillation frequencies were investigated and the 
so-obtained specific heats were extrapolated to zero 
frequency in order to obtain quasi-static values of the 
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specific heat itself. Specific heats of toluene and OTP 
[53 |. and PPhO [23[ were taken from the literature. 
Since the configurational entropy cannot be deter- 
mined experimentally, the entropy of the liquid in ex- 
cess to the corresponding crystal was used exploiting 
the proportionality between the two magnitudes Q- 
This only scales the value of C in Eq. (1) and a in 
Eq. (3). A linear form of the excess specific heat: 
AC P = C™ e " - C^ rystal , where a and b are material 
specific constants, was employed to obtain the excess 
entropy [2^|. This was obtained integrating the rela- 
tion: J^JAC p {T)/T]dT; where T K is the Kauzmann 
temperature, namely the temperature where the en- 
tropy of the liquid equals that of the crystal. Due 
to the lack of crystalline specific heat data for glass- 
forming polymers, the latter quantity was determined 
identifying it with the temperature where the dielec- 
tric relaxation time of the a process tends to diverge 
Hence, the knowledge of the parameters a and b 
allows fitting pure components dynamics data through 
the AG equation (dashed line in Fig. 3 for pure PVAc) 
to obtain C, log To and 7k. The values of the param- 
eters a and b together with T g , Tk, C and logr are 
summarized in table 1 for the pure polymers. Litera- 
ture values of the packing (l p ) and Kuhn (l^) length, 
relating the size of CRR to the self-concentration [16[ , 
are also included. In addition, previously obtained 
BDS and calorimetric data for PVAc in poly(cthylcnc 
oxide) [IH, for PoCIS in PVME and low molecular 
weight PS [H and PVME in toluene, PoCIS and PS 



with various molecular weight 15| were employed 



PVME PVAc 



PS 



PoCIS 



T g (K) 249 304 373 402 

T K (K) 203.5 256 322 340 

log to -13.0 -13.8 -12.4 -12.9 

Cpmol^KT 1 ) 21900 41450 26400 61800 

a(Jmor 1 K- 1 ) 68.3 105.3 99.0 98.4 

^Jmol^K" 1 ) -0.16 -0.21 -0.17 -0.17 

Z k (A) 13 17 18 17.3 

Z P (A) 2.7 3.7 3.9 3.9 

a(J 1/3 mor 1/3 K~ 1/3 ) 10.7 ±0.5 17 ± 0.5 19.5 ± 0.5 24 ± 2 



TABLE I: Relevant parameters for all components of all 
mixtures. Data for the C parameter of the AG relation, 
the parameters a and b of the excess specific heat and a 
are referred to a per mole of monomer basis. 

In what follows, we will present the results obtained 
for mixtures at high concentration of PVAc as a show- 
case and illustrate the methodology followed to eval- 
uate the parameter for this polymer. Fig. 2 shows 
the dielectric loss versus frequency for PVAc at vari- 
ous concentrations of toluene and fixed temperature. 
The increase in toluene concentration clearly provokes 
a speed-up of PVAc dynamics. Similar results were 
obtained for PVAc in PEO and BCDE. The average 
relaxation time of PVAc in mixtures with toluene as 
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FIG. 2: Permittivity loss vs. logarithm of the frequency 
for PVAc/toluene system at 308K and the following weight 
percentages of PVAc: 60% (filled circles), 63% (empty cir- 
cles), 71% (filled triangles) and 86% (empty triangles). 



well as that of pure PVAc as a function of the inverse 
temperature is displayed in Fig. 3. The continuous 
lines represent the best fits of the proposed model to 
experimental data. The quality of the fits is excellent. 
Analogously, the fits performed for the dynamics of 
PVAc in mixture with PEO and BCDE resulted in a 
successful description of data. The parameter a ob- 
tained from the best fits of our data is plotted in the 
inset of Fig. 3 as a function of the average effective 
concentration. Strikingly the extrapolation to pure 
PVAc leads to a single value of a independently of 
the environment surrounding PVAc chains. This can 
be regarded as the intrinsic a connecting PVAc con- 
figurational entropy with its size of CRR. The same 
procedure was applied to evaluate the parameter a 
for PS investigating its dynamics in OTP and PPhO 
and similarly to PVAc, a single a was obtained. Fur- 
thermore, a single value of a was obtained for PoCIS 
and PVME studying their dynamics in various envi- 
ronments. The so-obtained values of a are tabulated 
in table 1. 

The availability of the polymer specific parame- 
ter a allows determining via equation (3) the size 
of CRR and its temperature variation, starting from 
the knowledge of the excess entropy determined from 
calorimetric measurements. The so-obtained diame- 
ter of CRR is displayed in Fig. 4 as a function of 
temperature for the investigated polymers. Though 
the size of CRR clearly depends on temperature, its 
variation is rather smooth (« 20% for the diameter). 
However, we note that this can translate in a dramatic 
change of the number of particles involved in the coop- 
erative rearrangement. Furthermore, a clear relation 
with the rigidity of the chemical structure can be ob- 
served, being the size of CRR larger for the most rigid 
polymer, namely PoCIS, and smaller for the most flex- 
ible, namely PVME. Remarkably the addition of one 
bulky chlorine atom in the phenyl ring of PS results 
in a substantial increase in the size of CRR for the re- 
sulting polymer, namely PoCIS. These results provide 
a positive indication of the reliability of the approach 
followed. The dependence of the size of CRR on the 



4 




3.0 3.5 
1000/7" (K' 1 ) 

FIG. 3: logarithm of the relaxation time vs. temperature 
for PVAc segmental dynamics in PVAc/toluene systems 
with the following PVAc weight percentages: 60% (filled 
circles), 63% (empty circles), 71% (filled down-triangles), 
86% (empty down-triangles) and 95% (filled square) . Pure 
PVAc segmental dynamics data are also shown (empty 
up-triangles) . The solid lines are the fits of the model 
to PVAc segmental dynamics in toluene solutions and the 
dashed line is the fit through the AG equation to pure 
PVAc dynamics data. The inset of the figure represents 
the variation of the a parameter with the average effective 
concentration of PVAc in toluene (filled circles), BCDE 
(empty triangles) and PEO (empty circles). The solid lines 
are a linear fit of experimental data, which, extrapolating 
to pure PVAc, allow the evaluation of the a parameter for 
pure PVAc. 




FIG. 4: Size of CRR vs. the temperature normalized at 
the T g for all investigated polymers. 



rigidity of the polymer suggests a positive correlation 
between the size itself and the Kuhn segment, shown 
in table 1. The assessment of the size of CRR allows 
performing a rough estimation of the number of parti- 
cles involved in the cooperative rearrangement. This 
can be done considering the inter-chain distance for 
the investigated glass-forming polymers, which can be 
obtained from the static structure factor S{Q). Val- 
ues for this distance can be calculated from the avail- 
able measurements of the static structure factor of PS 
IS S3, PVAc [H and PVME [H|, being the resulting 
values: dps = 9A, dpvAc = 8.4A and dpvME = 6.3A. 
It is noteworthy that the diameters of CRR found at 
T g for these polymers are almost quantitatively twice 



the corresponding intermolecular distance, suggesting 
that the cooperative rearrangement only involves the 
first shell around a basic structural unit. This would 
imply that the number of particles involved in each 
CRR at T g is of the order of tens. The possible connec- 
tion between the size of CRR and the intermolecular 
distance opens interesting scenarios in the rational- 
ization of the length scale involved in the structural 
relaxation of glass-formers, which needs to be further 
explored. 

Finally, it is worth comparing the sizes of CRR ob- 
tained here with the AG approach with those previ- 
ously estimated employing different approaches. The 
size of CRR obtained for PVAc at T g , £ « 1.7nm 
agrees almost quantitatively with the characteristic 
length scale obtained by Berthier et al. [H by means of 
multipoint dynamical susceptibility (£ rj 2nm). The 
size of CRR obtained by us for PS is also in qualitative 
agreement with that estimated by Hempel et al. [29| 
in the framework of Donth approach Q (£ ~ 2.2nm vs. 
£ w 3nm) and in excellent agreement with that esti- 
mated by Sills et al. by means of nanoscale sliding 
friction (£ w 2.1nm). Our results are also compatible 
with the absence of finite size effects in confinement 
down to 5 nm (3l| . 

To summarize, we have taken advantage of the con- 
cept of self-concentration in combination with the AG 
theory of the glass transition to open a new route to 
obtain the size of CRR and its temperature depen- 
dence for several glass-forming polymers. This new 
approach proposed allows the unambiguous derivation 
of the polymer specific parameter a connecting the 
size of CRR to the configurational entropy. Once a 
is obtained for each polymer, the size of CRR in a 
wide temperature range can be determined from the 
knowledge of the configurational entropy of the pure 
polymer. The temperature variation of the diame- 
ter of CRR for a given glass-forming polymer results 
to be of only about 20% in the accessible tempera- 
ture range. Remarkably, the values of the diameter of 
CRR at T g for the investigated polymers closely cor- 
relate with both the Kuhn length and the inter-chain 
distance. 
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